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Abstract 

(N . 

■ We study the deformation quantisation (Moyal quantisation) of 

general constrained Hamiltonian systems. It is shown how second 

f*"^ \ class constraints can be turned into first class quantum constraints. 

This is illustrated by the O(N) non- linear a- model. Some new light 
is also shed on the Dirac bracket. Furthermore, it is shown how clas- 
i , sical constraints not in involution with the classical Hamiltonian, can 

be turned into quantum constraints in involution with respect to the 
Hamiltonian. Conditions on the existence of anomalies are also de- 

• , rived, and it is shown how some kinds of anomalies can be removed. 

rS [ The equations defining the set of physical states are also given. 

It turns out that the deformation quantisation of pure Yang-Mills the- 
ory is straightforward whereas gravity is anomalous. A formal solution 
to the Yang-Mills quantum constraints is found. In the ADM for- 
malism of gravity the anomaly is very complicated and the equations 
picking out physical states become infinite order functional differential 
equations, whereas the Ashtekar variables remedy both of these prob- 
lems - the anomaly becoming simply a central extension (Schwinger 
term) and the equations for physical states become finite order. 
We finally elaborate on the underlying geometrical structure and show 
the method to be compatible with BRST methods. 
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1 Introduction 



The problem of how to quantise a given classical system is one of the oldest 
in the history of quantum theory, its age being comparable to that of the 
measurement problem. This being so, many attempts have been made in the 
past of defining general quantisation schemes. Usually such schemes try to 
construct mappings qi,p* t— > fa, ft* such that 

{f(q,p),g(q,p)}pB ^[f(q,P),g(q,p)} (1) 

where {•, -}pb denotes the Poisson bracket, 

ff f \ / u df dg df dg 
{f(^9(q,p)hn- WiW - WWi (2) 

and where / denotes some kind of operator-valued function obtained from 
the classical observable / in some prescribed manner. The most elaborate 
of these schemes must be geometric quantisation, |]], which relies heavily on 
the symplectic geometry of classical phasespaces. 

Unfortunately, one can show that in general no quantisation procedure can 
exists for functions which are not at most quadratic polynomials in the basic 
variable, qi,pi, 0. This implies that the algebra will receive quantum cor- 
rections. In fortunate circumstances these are just central extensions. [] 
Another procedure, which has not been studied as much as geometric quan- 
tisation, is deformation quantisation. It is known that this always works, 
H (recently, Kontsevich has extended the proof of existence from symplec- 
tic to Poisson manifolds too, f27j). Denote the classical phasespace by T, 
this is then a symplectic manifold, and the set of observables, C°°(T), form 
a Poisson-Lie algebra (i.e., a Lie algebra under Poisson brackets), A. De- 
formation quantisation consists in replacing the algebra of observables with 
another, deformed one, An, and the Poisson brackets by a new, deformed 

1 One should note that BRST does not escape this problem either - given a classical 
BRST generator, the task is yet again to find an operator realisation. The problem is 
lightened a little bit, however, since BRST exploits the structure of the classical theory 
to a far greater extent than old fashioned canonical quantisation. Ultimately, though, it 
ends up with similar problems. 

Geometric quatisation avoids the no-go theorem by loosening another of Dirac's require- 
ments, namely that q,p be irreducibly represented. This can lead to problems with the 
classical limit, however. 
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bracket. More concretely, An = A<8>C((h)), i.e., the elements of An, are func- 
tions on phasespace which can also be expanded in a power series in h, fi" 1 .^ 
Thus elements of An can be written as 

f{q,p, H) = ... + rV-i^.p) + fo{q,p) + hh{q,p) + ... (3) 

where f n G A = C°°(r),n G Z. Here /o is the classical part of /. Similarly 
one writes the new bracket as 

[f,g] M = ih{f,g} PB + 0(h 2 ) (4) 

We will in general refer to the deformed brackets as Moyal brackets, even 
though that name is strictly speaking reserved for deformed brackets of flat 
phasespaces, T ~ M 2n . It has recently been proven by Tzanakis and Dimakis, 
0, that the Moyal bracket is essentially unique - the various possibilities 
corresponding to various operator ordering prescriptions. Consequently, de- 
formation quantisation does not suffer from the usual amibiguities of other 
quantisation schemes. Furthermore, Dereli and Vergin have proven that the 
various possibilities for the Moyal bracket du to different operator orderings, 



satisfy a Woo-covariance, p5| . The close relationship between H^-symmetry 
and the Wigner-Weyl-Moyal formalism has also been extensively studied by 
Gozzi and coworkers, ||26|| . 

Furthermore, we will usually restrict ourselves to the simpler case of only non- 
negative powers of h, i.e., work with An = A <8> C[[7i]]. Only when this turns 
out to be impossible will we deal with the more general case of A <E> C((h)). 
Systems with second class constraints or anomalies will in general require the 
presence of negative powers of h, whereas systems without such problems in 
general will not. 

The general philosophy of the paper will be that for any classical theory 
with problems (such as second class constraints, Hamiltonian not in involu- 
tion with the constraints, anomalies appearing upon a simple quantisation), 
a completely well-behaved quantum theory exists for which all of these prob- 
lems are absent. The problems only reappear when taking the classical limit, 
a procedure which has to be slightly modified. In other words, we have a 



2 Standard mathematical notation: C[x] denotes the set of complex polynomials in one 
variable x, C[[x]] the set of formal power series in x, C(x), the field of fractions of C[x] 
(i.e., C(x) = {f(x)/g(x) \f,g€ C[x],g(x) ^ 0}, and C((x)) the field effractions of C[[x}], 
i.e., C((a:)) ~ C[[x,x- X ]]. 
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"quantum smoothening out" of classical problems, analogous to what hap- 
pens in quantum chaos, where a chaotic classical Hamiltonian still leads 
to an integrable system at the quantum level, because the phasespace gets 
"smurred" by the Heisenberg uncertainty principle. I think this is a sensible 
procedure, as the quantum description is presumably the most fundamental 
one, and the classical description should at most be considered as a limiting 
case - and as such, we would expect problems to occur, not so much at the 
fundamental level (nature tends to favour simplicity) but arrising from tak- 
ing the limit. 

For T ~ R 2ra a compact expression for the Moyal bracket exists. Introduce 
the bidifferential operator A|] 

fAg = {f,g} PB (5) 

then 

\f,g] M = 2ifsm(±hA)g (6) 
see e.g. [0. This can be written out more explicitly as 

" (-1)^1 " / fl \ d^f d^g 

L^' y J M Z-, 4 n( 2n+ i)! Z-A ' y k J dq 2n ~ k+1 dp k dp 2n - k+1 dq k K ) 

with 

' n \ n! 



k J k\(n-k)\ 
denoting a binomial coefficient. 

The Moyal bracket can be written in terms of a so-called twisted product 

f*g = fei™*g = fg + 0(h) (8) 
the Moyal bracket is then the commutator with respect to this product 

[f,g]M = f * g - g * f (9) 

which brings out the analogy with the standard operator formulation of quan- 
tum theory. 



3 We will only in this paper deal with bosonic degrees of freedom. Fermionic variables 
can be treated by means of a Z2-grading, where the usual Poisson bracket is replaced by a 
graded analogue, as in the Batalin-Fradkin-Vilkovisky approach to BRST-quantisation. 
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In general we will write 



[f,9\M = ^2ih n u n (f,g) 



(10) 



with = {/,#} PB - 

Since one can always choose local Darboux coordinates, q%^ (with the stan- 
dard Poisson bracket) we can always at least locally make do with the Moyal 
bracket, without having to work with a more general deformation. Fur- 
thermore, for most systems of physical interest, such as all standard field 
theories, we can take the classical phasespace to be "flat" in the sense that 
it can be covered by a global coordinate patch of Darboux coordinates. The 
possible "curvature" of T is put into a set of constraints. Hence one may 
extend the original phasespace, which might not have a global Darboux co- 
ordinate patch, into a "flat" phasespace D, such that T is characterised by 
the vanishing of certain functions (f) a (q,p), i-e., 



This leads us naturally to study Hamiltonian systems with constraints. 



2 Hamiltonian Systems with Constraints 



Consider a classical Hamiltonian system with a phasespace T and a Hamil- 
tonian h(q,p). Suppose the transition from a Hamiltonian to a Lagrangian 
picture is singular, i.e., suppose not all the phasespace coordinates are inde- 
pendent. We then have a system with constraints, i.e., a family of functions 
0a(?,p) exists such that the true degrees of freedom are characterised by the 
vanishing of these functions. We will also assume these functions to be inde- 
pendent. 

We will, at first, assume that the set of constraints are first class, i.e., they 
satisfy a Poisson-Lie algebra 



where the structure coefficients, c ab c can be functions of the phasespace vari- 
ables. 

For now we will also assume that time evolution does not give rise to new 
constraints, i.e., the functions <p a are in involution with the Hamiltonian 



r = {(q,p) G T' | Va : <f> a (q,p) = 0} 



(11) 



{</>a, </>&}pB = C ab C (p c 



(12) 



{h,(f> a }pB 




4>b 



(13) 
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where again the structure coefficient can depend on the phasespace coor- 
dinates. We will return to the more general setting in a subsection later. 
Now we want to perform a deformation quantisation. We will thus replace <p a 
by functions $ a = <f) a + 0(h, Tr 1 ), the Hamiltonian h with H = h + 0(h, fir 1 ) 
and the Poisson brackets with Moyal ones, in such a manner that 

[<5>a,<5>b]M = ihc C ab <t> c (14) 

[H,* a ] u = ifiV^ b (15) 

For simplicity we will only consider Q a ,H to be formal powerseries in h and 
not fir 1 . Write 

oo oo 

$a = 0a + $> n $i n) H = h + Y,^H n (16) 

n=l n=l 

Furthermore we have already assumed that the structure coefficients are un- 
deformed. When the Moyal bracket algebra, (H-TED, become non-isomorphic 
to the classical Poisson bracket one we say we have an anomaly. This can 
happen if e.g. the structure constants get ^-corrections, the algebra become 
centrally extended or if the right hand sides of ( |I4T - |T5| ) become non-linear in 
the $ a . Such problems will be dealt with in a later section. 
At each order in h the Moyal algebra relations ( |T4} - |H)| ) read 

*? } ) = (17) 

n,k 

£uw*(ir B ,s«) = kM*° (is) 



n, k 



which provides us with a recursive scheme for finding the quantum constraints 
$ a and the quantum Hamiltonian H. Inserting the explicit form for u>i we 

get 



V-^" (-l)N-n-k + l f N _ n _ k 



4N-n-k(2N- 2n-2k+ 1)! V 1 

n,k 1=0 K 1 v 



X 



gN-n-k^( n ) QN-n-k^(k) 
dpN -n-k-lQ q l dq N-n-k-l dp l = '•<'< 



$(^19) 



Z-j 4 v-n-fc(2AT - in - 2k + 1)! \ I 



X 
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d N - n ~ k H n d N - n - k ® { a k) 
QpN-n-k-lQql QqN-n-k-lQpl 



V^i N \2Q) 



One should note that if <p ai h are all at most third order in q,p then u n , n > 3 
on these vanish identically and one can take simply $ a = (f) a , H = h. Many 
physical systems have precisely this structure. Most notably the Yang-Mills 
field theory. The study of membranes or some supersymmetrical theories, 
however, will in general need ^-corrections to the constraints and Hamilto- 
nian. 

It is well worth having a closer look at the recursive scheme above. The first 
equation, the one for the first quantum correction, reads 

c: b $( i ) = {0 a ,$i i) }p B + {$i i ),0 f) }p B (21) 

which gives a first order partial differential equation for , 



d d(p b d \ ( d(f) a d dcj) a d 



= (22) 



dq dp dp dq J \ dq dp dp dq / 
Introduce the operator 

Kb = c c ab + 5 c a {<j> b , -}p B - 6 c b {cf> a , -}p B (23) 
we can then write the equation for $c as 

= (24) 
The other equations in the hierarchy have a similar form, namely 

Kb®i n) = Jlb } (25) 

where J^b ls a source term only depending on $c = 4> c , &*c\ $c™ • 
Explicitly 

n m+1 

= ~ E E ^\ (26) 

m=l k=l 

which shows how the "source" is build up from the various m-order brackets, 

Note, by the way, that if ( 24]) has only the trivial solution $c = 0, then it 
follows from (25) that $£ = 0, n > 1 and hence $ a = a . This would of 
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course not be the case if we included negative powers of h in the expansion of 
$ c . In any case equations (|24| - f25| ) show that $c™' ) cannot be uniquely defined 
as we can always add a multiple of $£ to it (with an appropriate factor of 
h in front, naturally). 

Furthermore, the operator T> c ab cannot be identically zero, since the classical 
constraints satisfy 

K b <Pc = ~c ab C C (27) 

and are hence "eigenvectors" of V c ab . Equation (25) constitue a set of first 
order partial differential equations determining the quantities $i n ' ) . 
There is a cohomological element to this set of equations. Denote the algebra 
of constraints by g, q C A = C°°(T), and form the complex C n = A n Q of 
alternating n-linear expressions in the constraints (f) a . Introduce the operator 



Kb C : C 1 



C 2 by 



A abfc ■= {4>a, /b>PB - {06, fa}pB 

We extend this by linearity to a map C p — > C p+l , 



ab J c\...Op 



^2 si S n (°'){^(a)' U(b)a(c 2 )...a(c p )}pB 

a cyclic 



(28) 



(29) 



where the a is a cylic permutation of the letters a, b, C2, 
Then by construction 



.,c p , i.e., a e S, 



P+i- 



.4 
A 2 



C p 




c 



p+1 



Hence we have a cohomology theory 

Z p := Keri 

5 P := Imi 
H p := 



(30) 
(31) 



(32) 

(33) 
(34) 



By construction if 1 measures the obstruction to invertibility of A as an 
operator on functions, and hence to the uniqueness of $c • The cohomology 
cannot be entirely trivial since 



A 



ab 



2c„ 



(35) 
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i.e., the constraints themselves are a kind of eigenvector system for A. 
Now, the standard complex computing the cohomology of a Lie algebra is 
the Koszul complex which is precisely our C p complex, ||. The derivative 
operator, however, is given by 

sf(X\, X p+ i) = ^^( — l) 9+r /(Xi, X q , X r , X p+ i, [x q , X r }) (36) 
q<r 

where a hat over an x denotes that it is to be omitted and where [■, ■] denotes 
the appropriate Lie bracket (here the Poisson bracket). But the Poisson 
bracket is also a derivation, i.e., satisfies the Leibnitz rule 

{x, yz} PB = {x, y} PB z + y{x, z} PB (37) 

This is in fact one of the defining relations of a Poisson algebra, namely that 
it is a Lie algebra whose Lie bracket also satisfies the Leibnitz rule. From 
this one immediatly sees that A = s (possibly up to a sign), since 

{0, / 1 - p xi...Xp} PB = f 1 - p {(fi,xi...x p }pB 

= f 1 - p ^2{<fr,x q } PB x 1 ...Xq...x p 
g 

Consequently H* is nothing but the usual Lie algebra cohomology. 

In a somewhat compact, symbolic notation the condition (fI3| ) then implies 

that c$ and J differ by an exact term. Explicitly 

Cab*^ = {0a, <Mpb + exact form ( 38 ) 

etc. For Lie algebras with vanishing cohomology we then see that the higher 
order terms are completely determined by the Uk brackets on the lower 
order terms. This ensures existence of such quantum constraints. It further- 
more shows that equivalent quantum constraints only differ by exact terms, 
i.e., two set of quantum constraints are equivalent (give the same Moyal al- 
gebra) if and only if they are cohomologous (i.e., that they at each power of 
H only differ by an exact form). 

A closely related cohomology generator will appear again when we deal with 
second class constraints. 

This procedure of simply replacing Poisson brackets by Moyal ones has been 
used by Plebahski and coworkers in the study of self-dual gravity, large N 
limit of SU(N) theories and Nahm equations, f22fl . The only new thing here 



is the level of generality - instead of considering a particular theory with 
constraints, we here deal with a generic situation. 
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2.1 Some Particular Solutions 

It is quite natural to assume $ a ~ a (where ~ denotes weak equivalence 
in the sense of Dirac), or more concretely that $ a only depends on {q,p) 
through the classical constraints <f> a . Assuming this, the first equation in the 
recursive scheme reads 



d<5>£ ] d 

907 Cbc<pd ^ c 



(0) = c^-jrf - - cU d ^- (39) 



a solution of which is <&£^ = a" b 0a0& with a constant satisfying 

«? = «c a (40) 
2«fcL-2afcL = c>f (41) 

This set of equations cannot always be guaranteed to have a non-trivial 
solution. Abelian constraint algebras are obviously not a problem, so let us 
consider the simplest non- Abelian algebra 

{0i,<Mpb = 02 

in this case the structure coefficient is cf 2 = 1 and the equation for a"f 
reduces to 

(S a iS b 2 - S a2 S bl ) a 9 2 f = 2 (a 2 b f 5 al - a\ s 5 a2 - a 2 a f 5 bl + a]/ 6^ S g2 
which has the solution 

a 22 = arbitrary 
a\ s = a] 1 = 

12 1 22 

«i = -2 a 2 
whence the quantum constraints read 

1 ■ 



$1 = 0i + ft ^~2^1 + "102 J 02 
$2 = 02 + ft«2 2 02 

This then constitute an example with <3> a a , where the quantum con- 
straints depend on the phasespace variable only through the classical con- 
straints. One should note that such a construction has found use in general 
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relativity, where a new Abelian constraint replacing the Hamiltonian con- 



straint has been found, [14]. In that case, however, since the formal ma- 



nipulations were carried out at a classical level, Planck's constant did not 
appear as a deformation parameter. Nevertheless, it is interesting to note 
the similarity. We will return to this new gravitational constraint algebra in 
a later section. 

For the case of an su 2 algebra, {(fi a , (p }pB = 2ze afe c </> c , the conditions on the 
coefficients ot b „ c read 



a 



2/ x9 3/ eg . If eg 3/ eg 1 • of 

a 2 % ~ a 2 °2 + a l % ~ «1 Ol = ^ 3 

2/ r9 3/ eg If eg , 2/ eg 1 . of 



-a 1 / 51 + a 3 / 81 - off 5 9 2 + off 5{ 



gf 



which only has the trivial solution a b ^ = 0, illustrating that not all algebras 
can be treated in this way. The example of an O(N) non-linear a- model 
belongs to this category as will be shown later. 



2.2 The Physical States 

Now, this treatment of constrained systems have dealt with merely the kine- 
matics and not the dynamics proper, i.e., we have not studied the equations of 
motion. In the standard Dirac treatment of constrained Hamiltonian systems 
the first class constraints become operators upon quantisation, i.e., physical 
states have to satisfy 

0a|^) = O Va (42) 

where <f) a is some operator realisation of <p a . How is this modified in a defor- 
mation quantisation? 

A good starting point is to consider the deformation quantisation analogue 
of a state namely a Wigner function W^. We could then impose the 
constraints in the following algebraic way: 

$ a *W 9 = Va 

or more symmetrically (where [•, denotes the "anti-Moyal bracket", [/, g\\ [ 
f*g + g*f = 2fcos(±hA)g) 

[$a,W*]+=0 Va (43) 
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Notice, 

0=[$ o ,W*]+ =2$ a W* + ... (44) 

so to lowest order we have = away from the constraint surface. More 
precisely 

suppWf C f|Ker0 a (45) 
where is the h° component of W^. 

The physical Hilbert space, 7i p h y s = {|^) I Va : 4> a \^) = 0}, then gets 
replaced by the space 

C phys :={W\ Va:[$ a ,W] + = 0} (46) 

This also has the advantage of being analogous to the BRST-condition, [Cl, A] = 
where [•, •] is a graded commutator, (l is the BRST-charge, (l = 7] a cj) a + ... and 
A is any observable (the observable corresponding to a state \^>) is of course the 
projection operator A = |^)(^|). Thus we will be using ( f43|) as defining the de- 
formation quantisation analogue of the Dirac condition. 

A few comments are in order. First, the replacement of Poisson brackets by Moyal 
ones implies that the corresponding "gauge" transformations acquire quantum 
modifications. If the classical constraints are cp a , then they generate (infinitesi- 
mal) "gauge" transformations 5 w f := {uj a (p a , /}pb> the corresponding quantum 
version is 

8 W F := [uj a <5> a , F] M = ih{u a 4> a , F} PB + other terms (47) 

which a priori differs from the classical expression. The discrepancy between the 
classical and the quantum "gauge" transformations show up in higher order deriva- 
tives, which seems to suggest that the quantum transformations are "larger", i.e., 
slightly less local than their classical counterparts. 

Another point to check is whether the space of physical quantities is invariant un- 
der such transformations. Consider thus an element A satisfying [$ a , A]^ = 0, Va, 
when then wants to prove that a "gauge" transformation does not take us away 
from this subspace, i.e., 5 u [$> a , A]\ t = 0, Va. We get 

(L([$ a ,A]+) = [<5 W * ,A]+ +[* 0> 5 W A]+ 

= [io b [^ b ^ c }M,A]+ +[$ a ,[w b *b,A] M ]i { (48) 

By noting that the physicality condition implies $ a * A = —A * <£ a , Va, we can 
then rewrite this as 

A]+) = J ([[$ 6 , $ a ] M , A}+ + [[* BJ A]+) = (49) 
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Hence the condition = [<& 0)J 4]]{^ is a reasonable quantum analogue of 4> a = 0, 
as we had anticipated, since it fulfills the only two requirements on can make a 
priori, namely that it has the correct classical limit (the correspondence principle) 
and that it is invariant under (quantum) gauge transformations. Notice that this 
even holds in the case where the structure coefficients depend upon the phasespace 
variables, as e.g. in gravity. 

For completeness we recall how Wigner functions are formed. In standard quantum 
mechanics in d dimensions one defines [|| 

W^(q,p) = J ^{x + ~y) ® i>(x - iy)e- i *P J??U (50) 

where ip is some wave function. This definition is appropriate in flat phasespaces 
and in the absence of external gauge fields. The case of gauge fields, however, can 
be treated by defining ip(x ± ^y) in a covariant manner using parallell transport, 
||. This can also be extended to quantum theory in curved spacetimes, [j/J. More 
generally, one defines by means of a Wigner- Weyl-Moyal (WWM) map as 

T^ = Tr(I%,pM<V>|) (51) 

or for mixed states in terms of a density matrix p as W p = Tr(U(q,p)p(q,p)), 
where II encodes all of the deformation information, ||], i.e., an operator A gets 
mapped to a phasespace function Aw (its Weyl symbol) 

A w = Tr(ni) (52) 

and where the deformed product is given by 

A w * B w = Tr(ILiB) = (AB) W (53) 

and so on. It has been proven that such a WWM-map II exists for a very large 
range of algebraic structures, @. 

One of the main features of the WWM formalism is that states and observables 
are treated on an equal footing, in fact is given by a "projection operator" in An, 
i.e., an *-idempotent P, P * P = P. Thus C p h ys C An- 

This can be clarified a bit by noting that the set of possible density matrices p, 
and hence of possible Wigner functions, is the (closed) convex hull of the set of 
*-idempotents (projections), 

C = co{P eA h \P*P = P} (54) 

and consequently, C p h ys C C. 

As a final comment worth making, we could note the possibility of letting the 
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quantum constraints be the classical ones but with the usual product of phasespace 
variables replaced by the twisted one in a symmetric manner, i.e., replace qp 
by ^(q * p + p * q). But when these are Darboux coordinates, i.e., when they 
are canonically conjugate, then \{q* p + p * q) and qp = pq are identical, since 
q*p = qp+ih,p*q = qp—ih. It is only when we are working with a phasespace which 
is not covered by a global Darboux coordinate patch that this possibility have any 
relevance. Since the non-flatness of the phasespace can always be absorbed into an 
appropriate set of constraints on a larger, flat phasespace, this possibility is only 
of academic interest. 

3 Overcoming the Various Problems 

At the classical level, a number of problematic features can be present. The con- 
straints may not all be first class or the constraints might not be in involution 
with the Hamiltonian, in both cases time evolution will take one away from the 
original constraint surface. We will show how to overcome these problems within 
the formalism of deformation quantisation. More precisely, we want to show that 
provided one is willing to allow negative powers of h (which is standard in the math- 
ematical treatment of deformation quantisation) , one can remove at the quantum 
level most if not all of the problems present at the classical level. It is this which 
is the main result of this paper. 

Problems may also arrise after a naive attempt at quantisation, such problems will 
be referred to as anomalies, even though they may not be just the familiar cases of 
central extensions of the algebra (a.k.a. Schwinger terms in current algebra) but 
can encompass a wider range of problems. These too can be overcome at least in 
some cases as will be shown. 

3.1 Second Class Constraints 

In Dirac's terminology, second class constraints are constraints i/j a which do not 
have weakly vanishing Poisson brackets, i.e., 

{ipa, A}pb = C ab 96 det C + (55) 

where being weakly zero means vanishing when the constraints have been taken 
into account (i.e., vanishing on the constraint surface). Dirac originally dealt with 
second class constraints by defining a new bracket, the Dirac bracket, 

{/,#}db = {/,#}pb - {/,^a}pB(C _1 ) a6 {5,^ 6 }p B (56) 
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which then have the wonderful property 



Wa,^}DB = (57) 

i.e., the bracket of two second clas constraints become strongly zero. 
There are a few problems with the Dirac bracket, however, most notably that it 
is usually very hard to compute explicitly. As far as deformation quantisation 
is concerned though, one can simply replace Poisson brackets by Dirac ones in 
the definition of the Moyal bracket to get a modified Moyal bracket. We will not 
pursue that line of inquiery any further here, since this approach will suffer even 
more from the calculational difficulties of the Dirac bracket. Instead we will ask 
two other questions: (1) is it possible to define "first class" quantum constraints 
^a(K) with ^ a = ... + ip a + + ••• (we will have to include negative powers of 

h in order to be able to cancel the C ab from the zero'th order term), and (2) can 
the Dirac bracket itself be seen as a deformation of the Poisson bracket? 
Let us first look for quantum constraints ^ a satisfying 

[*o, *&]m = ifa&tf c (58) 

where we explicitly want to be able to have e c ab = too. Since we now have to 
include negative powers of the deformation parameter h, we write 

00 C ) c 

¥ a = X>*j: B) < 6 = 5> n e a6 (59) 

raeZ n=l 

and compute the Moyal bracket. The ftP part of this becomes 

oo c oo 

E ^ n ~ 1} = EM n) '*i" n) }pB+EE^(^ n) '*i" fe_n) ) (so) 

n=l nel k=l rteZ 

which is then one of a family of equations the new constraints have to satisfy. 

(n) 

It will often be enough to assume \l/c = 0, n = —2, —3, —4, i.e., that only one 
negative power of the deformation parameter occurs. In this instance the set of 
equations read 

= W^V^Ipb, (61) 



= {^V^PB + WvV'i >B, (62) 

{ e C ab ^~ l) -C ab = a; 3 (^- 1) ,^" 1) ) + M- 1) ,tf ) }PB + {V'i 1) ,4" 1) }PB, 

(63) 
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and so on. One would often like to take e c ab = but this might not always be 
possible, and not all choices of structure coefficients may be allowed. There will, 
however, be quite a lot of freedom in the particular choice of algebra in the general 
case. 

As was the case for first class constraints, we can introduce a cohomology theory 
for this set of equations. Define 

B C a lfc,., n = £ S[ ^M~(al /^)--(cn)}PB (64) 

then B 2 = 0. Notice that this is the same cohomology construction as for first class 
constraints, except that we now use the tr l components and not the semiclassical 
ones for the definition of the differential. With this the conditions read 

K b 4~ 1] = (65) 
K^c = (66) 

K h 4 l) = i e C ab 4- 1) +C ab -^i- 1 \4- 1) ) (67) 

So, contrary to the first class case, both the the h~ l and the h° terms are closed. 
Let 

Z p = KevB B p = lmB H P = Z P /B P (68) 

then 

^^aez 1 (69) 

but 

^ez 1 !^ 1 ^ 1 ) (70) 

which shows how the subsequent terms are build up from cohomological ingredient 
at lower levels, i.e., that the h-teim is a closed form plus the third order bracket, uj^, 
of two closed forms. Hence, as for the case of first class constraints, the solutions 
are described by a cohomology. For second class constraints, however, the algebra 
of constraints is somewhat different. Instead of simply the algebra g C A = C°°(T) 
of classical constraints, ip a , we have an extended algebra gC^ = C°°(T), which is 
an extension of g by an Abclian algebra go spanned by the h" 1 components. The 
algebra go defines a cohomology, H p , whose differential operator, B, extends to 
all of the extended algebra g, making the situation somewhat more difficult. The 
extension g need not be a central extension, since the condition B^ b tfj c = does 

not imply {V'ajV^ ^}pb = 0. It is clear, however, that central extensions would 
be a good starting point when one wants to construct g. This once more brings us 
back to the Lie algebra cohomology of g, since the second cohomology class thereof 
describes the various possible central extensions. We will not elaborate more on 
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the cohomological element here. 

It is important to emphasise that the desire to remove second class constraints 
forces the pressence of negative powers of h upon us. But provided one is willing 
to pay this small price, one can treat second and first class constraints in the same 
manner at the quantum level. This of course implies that one has the redefine 
the classical limit slightly, a point to which we will return towards the end of this 
paper. 

In the following subsection we will consider a particular example of a physical 
system with second class constraints and show how to use deformation quantisation 
explicitly. But first we will return to the Dirac bracket construction again. 



The Dirac bracket, (56), can be seen as a deformation of the classical Poisson 
bracket in its own right. The deformation parameter is in this case a matrix 
valued function, namely C ab itself. This is most easily seen in the following simple 
toy model where we have an algebra of first class constraints (j) a and a pair of 
second class ones, ip,ip, satisfying the generic algebra 

{^aAb]pB = c c ab (f) c 

{4>a,4>}pB = C b a (f) b 
{4>a,$}pB = C b a (f) b 

{ip,ip}pB = k 
The Dirac bracket of the first class constraints then reads 

{4>a, 0fe}DB = C C ab (j) c - k~ X C C a ^ c ^ d 

One can either see this as defining a "quadratic algebra" or as a deformation of 
the Poisson bracket. This latter point is particularly clear when one considers two 
arbitrary functions /, g and compute their Dirac bracket 

{/, s}db = {/, £/}pb - k^itp, /}pb{^, g}pB ■= {/, g}pB - fc~M/> g) 

Contrary to the Moyal bracket deformation of the classical Poisson brackets, the 
Dirac bracket, however, does not involve higher order derivatives of the observables. 
It is a deformation all the same; the limit k — > oo corresponding to the classical 
Poisson bracket. Contrary to the Moyal bracket, the higher oder brackets do not 
involve higher order derivatives. Deformations such that the n'th order bracket 
contains only n'th order derivatives is known as a deformation of Vey type. Hence, 
the Moyal bracket is of Vey type whereas the Dirac bracket is not. 
The Dirac bracket amounts to replacing the Poisson bracket by a new first order 
bracket 

{/, gh = ^(C)-- - VAC)-- - MC)-- + a PP (C)^ Yp (71) 



17 



where 



a 



(C) = 0(C) a qq = 0(C) (72) 

P pq (C) = -1 + 0(C) P qp = +1 + 0(C) (73) 

where C is allowed to be a matrix valued function of the phasespace variables. We 
will use the name secondary deformations for such deformations. A priori, such 
secondary brackets will not define Lie algebras; one can only ensure this to 0(C). 



Before we move on to consider some examples and to study anomalies, it is worth 
recalling that other ways of removing second class constraints exist. These are 



inspired by the BRST approach to quantisation, [18|. One enlarges the phasespace 
by including Grassmann odd variables, rj, V, and then replaces the constraints tpi 
by VPj = ipi + X^n=i X-'Xj where x can be either r\ or V . The extra Grassmann 

n 

variables (the ghosts) here correspond to our deformation parameter H, and is 
just one indication of a connection between BRST and deformation quantisation. 
We will comment a bit more on this connection in the section on the underlying 
geometrical structure. 



3.1.1 Example: O(N) Non-Linear a-Model 

Consider an 0(iV)-valued field 4>, and take the Lagrangean to be 

£ = X -d^ A d^A + \\(4> A 4>a - 1) (74) 

where ^ is a (flat) spacetime index and A is the 0(iV)-index. In the Hamiltonian 
formalism of this O(N) non-linear cr-model, we have the following constraints 

V>i = 0A0 A -1 (75) 

i>2 = TT^A (76) 

where it a is the canonically conjugate momentum of 4>a- The constraints are 
second class and satisfy 

(ViOr), Mx')}pb = 25(x - x')4>a4> A (77) 

4 This is somewhat similar to the situation in equivariant cohomology - if one attempts 
to use the secondary bracket {•, -}2 to define a cohomology theory, i.e., a derivation s, one 
will get s 2 = 0(C) in the generic case, and one would not get a standard cohomology 
theory. 
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which is clearly non-zero on the constraint surface. 

We now want to find quantum modified constraints = 1,2 with vanishing 
Moyal brackets and with 

= nr V{ _1) + & (78) 

This leads to the following set of equations (suppressing the spacetime variables 
x, x' and the resulting Dirac 5- functions) 

{^ _1) ,^ _1) }pb = (79) 

{Y^Y^pb + IY^V^Ipb = (80) 

= {i>i,i> j }pB = 2e ij <t> A <t> A (81) 

uj n {4~ 1] = n>5 (82) 

We will make the following Ansatz for the quantum part of the constraints 

4 _1) =^i/W) 4 _1) =^M<M (83) 
This yield the following set of coupled equations 

= 2^ 2 /5 + V'i#-|| + ^iV'2||-^ + 2/^-|^ (84) 
= 2<t ) 2 (f + g) + 2^.^- + ^-^ (85) 

= i 25AB 4l- + 2 ^ fifL + ^ " + c y clic l x 

V 9</>c d(f>BO(f>c 3 d(f> A d(f>Bd(f>c ) 



ftrW + 3^ • "^w^ + cyclic J " 2<p (86) 

where cyclic denotes a sum over cyclic permutations of A, B, C. 

We can ensure w„(fi,f2) = 0, n > 5 if we take g to be only quadratic in the 

momenta, i.e., 

g(<f>, tt) = + /? A 7T B + 7 (87) 

where a A B,(3 A ,7 can depend on 0. Notice, furthermore, that the (^-derivatives of 
g do not appear at all in this set of quations, hence there is a lot of freedom in 
choosing the ^-dependency of g. We will use this to assume a, [3 to be independent 

of (j). 

By combining the first two equations we get 

24> 2 fg + 2f(<f>-ir)<f> A (a AB Tr B + (3 A ) = - f (<t? - 1)4> ■ d f (88) 
(g -f)<j>. df + ^-7r)(a AB 7r B + (3 A )d A f = (89) 
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where we have written df for K . Collecting powers of the momentum we get 

a A BK A K B 4> c dcf = -a B cK A K B 4> c d A f (90) 
pA<f> B n A d B f = -pA<p B 7r B d A f (91) 
( 7 -/)^ A / = (92) 

To get an understanding of these equations we will consider first the case of N = 1. 
The third condition then implies 7 = /. Using this we arrive at 

2/' + 2#" + (0 2 -l)/ w = 20/a (93) 

The form of this equation suggests an expansion on Hermite polynomials of h = f. 
Thus, write 

00 

h = Y J C n H n (<P) (94) 

n=0 

Inserting this in (|93|) we get 



^2 [cn + 20(n + l)c n+1 + 2(n + l)(n + 2)(</> 2 - l)c n+2 ] # n = — #x (95) 

n=0 

where we have used the explicit form for H\ (i.e., H\{x) = 2x) and the standard 
formula 

—H n (x) = 2nH n -i(x) 
ax 

aii' 

ity of the Hermite polynomials 



Now, multiplying by H±e ^ and integrating over cj), we get from the orthonormal- 



2^F [a + 12c 3 - 12c 3 + 480c 5 + 72c 3 ] = — (96) 

a 



From this we get 



h = cxH x + c 3 H 3 + C5F5 (97) 
= (2ci-12c 3 + 12Oc5)</> + (8c3-16Oc 5 )^ 3 + 32c50 5 (98) 

where the coefficients ci , c 3 , C5 are restricted by (^) . There are two undetermiend 
coefficients since the equation defining h is of second order. The function / is then 
the primitive of h, i.e., 

f(4>) = c + (ci - 6c 3 + 6Oc 5 )0 2 + (2c 3 - 40c 5 )^ 4 + yc 5 </) 6 (99) 
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where Co is an undetermined constant. Note, by the way, that, on the constraint 
surface, / = const. Hence our new quantum constraints are 

*! = ^(i + zt 1 /^)) (ioo) 

*2 = ^(l + /T 1 (a7r 2 + /37r + /($)) (101) 

We have now constructed a set of quantum constraints ^i, ^2 satisfying an abelian 
algebra under the Moyal bracket, 

[*i,* 2 ]m = (102) 

and the second class constraints have consequently been "lifted" to an abelian, but 
deformed, algebra. This illustrates the general procedure of deformation quantisa- 
tion of Hamiltonian systems with second class constraints. For the remaining part 
of this paper, we will implicitly assume that second class constraints have been 
dealt with in this manner, and will thus concentrate on first class constraints. 
The case N > 1 is treated in an analogous manner, but the solution will here be 
in terms of many-variable Hermite polynomials, -ffm,...,n JV (0i) ■••> <A?v) which satisfy 
similar recursion relations. Unfortunately, the explicit expressions quickly become 
highly complicated and messy and we will not give them here. It is, however, 
fairly straighforward to do so - it merely corresponds to finding solutions of the 
many-dimensional harmonic oscillator Schrodinegr equation. In a multi index no- 
tation, moreover, the equations essentialy reduce to the N = 1 case (with a few 
subtleties). 

3.2 Constraints not in Involution with the Hamilto- 
nian 

Another problem which may occur is the failure of the constraints to be in invo- 
lution with the Hamiltonian, i.e., 

{Mo}pb^K 6 & (103) 

In this instance one cannot simply impose the constraints for some initial values 
of the phasespace variables, but have to impose them at each time step, since time 
evolution does not preserve the constraints. 

Explicitly, we will assume the constraints to be first class but satisfying 

{h, <Mpb = Vfa + X a (104) 
where Xa 7^ 0. The Poisson bracket of Xa with the constraints will be denoted by 

{4>a,Xb}pB = Kb (105) 
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with A a f, 56 e c ab Xc for some function e c ab . This implies that Xa cannot be seen as 
merely a new first class constraint. 

The task is to find new quantum constraints and Hamiltonian such that 

[H, $Jm = ihV^ b (106) 

with = Va+0(h). We will again assume to be independent of the phasespace 
coordinates. 

Inserting an expansion in h, 

H = ^h n H n 

n 

$ a = 

n 

(n) 

n 

we get for the h° part 

(n) 

E <w 2 „" 2n+1 (tf fc , = E ^ ( 107 ) 

fc,Z,n n 

which clearly shows that negative powers of h are needed - otherwise we would 

(o) 

get the condition V£(f>b + Xa =V£ fib which is impossible. It does, however, make 
sense to assume that H and V% have only non-negative powers of the deformation 
parameter, whereas the only negative powers of h appear in 3> a . We can even 
assume that only powers h n ,n > —1 occur. 

With this assumption, the condition coming from the h~ l part reads 

(o) 

^{hMa^) (108) 

i.e., the corresponding term in the quantum constraints are in involution with the 
Hamiltonian, with the structure coefficients being given by the involutive part of 

(o) 

the classical relation, i.e., V£= V^. 

The "classical" part, i.e., the H° contribution, reads 

(o) (i) 

vfy b + xa + $ (- 1 )) =v h a 4> b + v a b *< x) (109) 
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(0) 

which simplifies upon putting V b = V£ as mentioned above. 

Since there are too many indeterminates in this problem (the H, <3? a and V b ), we 
will make a further simplifying Ansatz, namely H = h, i.e., that the Hamiltonian 
receives no quantum corrections at this level - all quantum modifications will be 
put in the constraints and the structure coefficients. With this condition we arrive 
at 

(i) 

V b a $ = X a (HO) 



Together with the involutive requirement, ( 108 ), this leads to a differential equation 

(i) 

for V£, namely 

(i) (i) 
{h,(V- 1 tx b }PB = V^V- 1 Y b Xc (HI) 

(i) 

Letting v b a =V b and using a matrix notation, this can be written as 

{h, v-^pbx + v- l {h, Xb }pB = Vv^x (112) 
which we can also write as 

vVv~ l x ~X = v{h, v~ 1 }pbX = ~{h, v}pbv' 1 x = -vv^x (113) 
In the simple case where we have only one constraint this equation becomes 

v d ^ d i 

V -dt lnx = -dt lnV 

which has the solution 

v(t) = c e- vt X 

where we have used that V is assumed independent of the phasespace coordinates 
in order to carry out the time integration. 

Hence constraints which at the classical level aren't in involution with the Hamilto- 
nian can, by deformation quantisation, be turned into quantum constraints which 
are. The price being that the constraints have a singular classical limit H — > 0. 
This price turns out to be the standard fare when dealing with classical con- 
straints with problems, or when one wants to eliminate problems arrising from a 
naive quantisation. 



3.3 Anomalies and their Lifting 

Upon deformation, the constraint algebra can be modified in a number of different 
ways. We will refer to such modification in general as anomalies. We have a 
natural hierarchy of such anomalies: 
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• Anomaly of zero'th order: [cf) a , (j) b ] M / ih{(f>a, 4>b}pB = ^c^ c . 

• Anomaly of first order: [$ a , <&b\M / ^ c a6^c 

The anomalies of zero'th order will in general be removable by deforming the 
constraints. One particularly important situation is when 

u n (4>a, 4>b) = k$ = const. (114) 
the Moyal bracket of two undeformed constraints is then 

[(pa, (pb]M = itlCab&c + Kb (115) 

with k ab = J2 n k^h 2n+1 . Hence we have a central extension of the original con- 
straint algebra. This is the way an anomaly very often shows up. We want to 
find deformed constraints, i.e., quantum constraints, <I> a such that the anomaly is 
"lifted", [& a ,&b]M = ihCab^c ~ wnere we allow for a change in algebra for h / 0. 
One should note that this procedure is analogous to the way we "lifted" second 
class constraints and the way in which we handled constraints not in involution 
with the Hamiltonian; the procedure thus shows the great flexibility and power of 
deformation quantisation. 

( n ) ( _ k) 
It will in general be sufficient to assume c c ah = for n < and ® a ' = for k > 2. 

(n) 

Furthermore, suppose k ab vanishes after a certain stage N, it is then natural to 

H (n ) 

choose c c ah = <& a = for n > N too. Moreover, this will be the typical situation as 
the constraints usually are polynomial and hence have vanishing n-order brackets 
for n sufficiently large. 
The general conditions become 

<4<^ -1) = U^aM^+Uli^Ab) (H6) 

2b + cl <P C + 2b ^ = "l(<t>a, $ I 1} ) + "l^, 4>b) + LO^,^) + 

Wl^j" 1 ),^) +U J3 ( ( f )a ,^ b ' 1) ) +L0 3 (^- 1 \4>b) 

(118) 

C C ab <S>^ + C C ab C C ab <f>c+ C C ab ^ = ^(d>( 2 ), fo) + U^a, ^) + U^,^) + 

u^-V, * ® ) + , ^ ) + A$ (119) 

and so on. 

We quickly recognize a cohomological element to these equations. In terms of the 
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aforementioned cohomology operators A^ b , B^ b , we can write the set of equations 
on the form 



(0) 

c 



c 



ab 



^ = i^-D = B c ab <f> c (120) 
( cL = (121) 

( cl C + cl = + ^*< 2 > + c 3 (^, *r i} ) + ^i-v^ b ) 

(122) 

etc., stating that the anomaly prevents <3?i ^ from being A-closed, or equivalently, 
4> a from being S-closed. In this way, an anomaly is recognized as an obstruction, 
and hence as corresponding to a non-trivial cohomology class. 
Let us again for simplicity restrict ourselves to the simplest possible non-Abelian 
algebra c\ 2 = 1, a,b,c = 1,2. We will then write simply k for k\2- We can 
furthermore assume tU2n+i(<Pai <Pb) = f° r n > 2, i.e., N = 2 (a general value for N 
can be treated similarly but the notation quickly becomes cluttered). This implies 

(n) 
L ab~ 

(0) 

c ab = c ab- The f un se l °f conditions read 



we can take c c b = ^> a n ^ = 0, n > 2 too. Moreover, it is natural in any case to pick 



^ = Wl ($S- 1) ,0 2 ) + a; 1 (0 1 ,^- 2) ) (123) 
cl^ l) = u l {<b[- 1 \^ ) )+uj 1 {^ ) M 2 ~ 1) ) (124) 

w 3 (*i _1) ,^2) (125) 

ci&P = Ul (*P,*P) + k (126) 

ustfu^ + uspQ-^th) = (127) 
It turns out, that a solution can be found if we make the Ansatze 

^)=a b J b + p a ^=7^ (128) 

where a, (3, 7 are constants. The appearance of a (3 is needed because of the fourth 
relation above would otherwise yield k oc (p which is explicitly assumed not to be 
the case. Furthermore, the first relation shows that $i ^ can only be proportional 
to <p a and not be linearly dependent upon it. 
By inserting the Ansatze in the relations above, one gets 

72 = 71=0 
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(1) 

12 12 12 

C12 7i = « 2 7i + «i72 
(i) 

Tra = a|+c? 2 
(i) 

a\+c\ 2 = 
ft = 

(i) 2 

det a = c 12 a c 
(1) i 

c\ 2 a\ = 
(i) 

c 12 Pc = k 

From which we get the structure coefficients 



and the linear part, (3, to be 

1 



ft = 1 72 — rr ^ = 72* (130) 



Which shows, that it is the central extension, k, which necessitates the linear part 
(3. Finally we get for the matrices a, 7 that they must be of the form 

«=( ai ~ ai ) 7=f ' 71 ) ( 131 ) 

^ q 2 -a x ) V 72 / 

subject to 

- 2a 2 7i = 0172 (132) 

but otherwise they can be chosen at random (but non-zero, otherwise 7 = or 
a = which is clearly not what we want). 
Explicitly, the new quantum constraints read 

$1 = ^ 1 7i02 + ^i + toi(</'i-02) + ^(72-ar 1 ) /c ( 133 ) 

$2 = ^ _1 72</ , 2 + 4>2 + fr(a2<fil - OL\(\>2 + 72&) (134) 

and the Moyal algebra is 

[$i,$ 2 ]m = (1 - a 2 h)$ 1 + CKift$ 2 (135) 
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which is still isomorphich to the classical, undeformed Poisson algebra, since that 
was the unique non-Abelian algebra with two generators. In general the Moyal 
and the Poisson algebra needn't be isomorphic, however, the latter only being a 
limit (or contraction), H — ► 0, of the former. 

The important thing to learn from this example is that the quantum constraints 
acquire quantum modifications of negative order in h as well as of positive orders, 
and that, furthermore, the structure coefficients too get quantum modified, but 
only with positive powers of h. Note, by the way, that it would be inconsistent to 
choose c^ b = c^ b , hence anomalies of the zero'th order can be lifted, provided one 
is willing to modify the structure coefficients. 

In the case where k ab is not a constant, we will have to let a, (3, 7 too depend on 
the phasespace variables. In this way, the conditions on these coefficient become 
differential equations, just like they did for the case of constraints not in involution 
with the Hamiltonian, and we will not comment further on that here. 
Anomalies of first order can result from a bad choice of quantum constraints, but 
if that is not the case, they cannot in general be lifted. It is difficult to think of 
any examples of such systems; they would correspond to a set-up in which the only 
consistent quantisations were drastically different from their classical counterparts. 
If the Moyal algebra becomes simply a non-linear algebra, i.e., if for instance 

[$a, &b]M = iftc ab c <& c + ih 2 d ab cd $ c <$> d 

then we can lift this anomaly by extending the constraint algebra to include Q ab := 
M> a <Iv Typically, this would lead to an infinite dimensional constraint algebra. 
Hamachi, |20[j , has recently considered another kind of anomaly "smoothening" . 
Again, it turns out that anomalous contributions can be removed at the Moyal 
algebra level. However, Hamachi restricts himself to systems with constraints 
only depending on momentum. The set-up presented here is much more general, 
although our results are not as rigorous as his. Just like in our case, the anomaly 
reappears as a troublesome H — ► limit. 

4 Examples: Yang- Mills and Gravity 

Possibly the two most important kind of theories in modern theoretical physics 
are Yang- Mills theories and General Relativity together describing all known forces 
in the universe. It is consequently of interest to devote some time to the study 
of their deformation quantisation. Especially when one considers the important 
differences in the structure of their respective constraint algebras. Where the 
structure coefficients of Yang-Mills theory are independent of the fields (thereby 
forming a true Lie algebra), the corresponding coefficients of the gravitational case 
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do depend on the fields (the metric to be precise). Furthermore, at least two set 
of constraint algebras exist for general relativity, the ADM ones and the Ashtekar 
variables, the latter formally resembling Yang- Mills theory quite a bit. Recently, 
a third algebra of constraints for gravity has been proposed, where the algebra is 
of the form C°°(£) x Diff(S) where E is the Cauchy hypersurface and x denotes 
semidirect product. We will return to these after having dealt with Yang-Mills 
theory. 

4.1 Yang-Mills Theory 

For a Yang-Mills theory with some gauge algebra g, the constraints are 

Qa = A< (136) 

where ir l a is the momentum conjugate to A?, Di denotes the gauge covariant deriva- 
tive and i is a space index (which we will take to run from one to three) and a is 
a Lie algebra index. The constraint algebra is 

{G a (x), Q b (x')} PB = c c ab 5(x - x')G c (x) (137) 

where c c ah are the structure coefficients of q. 

Since the constraints are quadratic in the variables (they have the form Q ~ 
8tt + ttA), we can take Q a to be the quantum constraints too, i.e., 

[Ga{x),g b (x')] M = ihc c ab 5{x - x')g c {x) (138) 
The last bracket is the one with the Hamiltonian 

n = \« - \f*F^ (139) 

Although TL is fourth order in the fields (it has an A 4 term) the constraints are 
merely quadratic and hence co n (H,G a ) = for n > 1. Thus we can take TL to 
be also the quantum Hamiltonian, and we get no anomalies (i.e., no anomalies 
from the deformation quantisation procedure, that is, other types of anomalies, 
especially global ones, are still possible). 

The constraint equations become finite order functional differential equations 

= [Qa, W] + M = 2(Dy a )W + lihH{cl b ^^ (140) 

where we have used [/, g]^ = 2/ cos(^A)g,(l(|. We have also suppressed gauge 
indices on the Wigner function (it has two - it taking values in g (g> g*, || [?]]). 
Formally, a solution can be found of the form 
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where g l i = <7 y (7r) is a quadratic function with values in the Lie algebra, g % i = 
d bc TT b T^ c where d bc is antisymmetric in be and satisfies 

c c ab d ba = 8ih- 2 (142) 

where indices are raised and lowered by the Kronecker delta. Since we already 
have one invariant antisymmetric tensor on q, namely the structure coefficient, 
c c ab , it is natural to put 

d c ab = aCc ab (143) 

which leads to 

a = - * (144) 

where k = 5 ab K a b is the trace of the Cartan-Killing metric, = c c ad c d bc . For 
= su n for instance, we have k = —2, leading to 



W[A,ir] = exp (-4ih- 2 J c bc F^iiTidx 



For SU2, teh Wigner function then involves the dual of the field strength tensor, 
whereas in other cases it involves some kind of generalised dual. In any case, the 
Wigner function is a kind of "bi-Gaussian" , i.e., a Gaussian in either A or ir with 
the coeffients depending on the conjugate variable. 

It is interesting to note that the present solution cannot be the Wigner function 
of a pure state. This is seen as follows. First, the Wigner function has the form 

W(q,p) = e iaq+bq2)p2 

hence the Fourier transform p — > y looks like 



W(q,y) = J * e 
V aq + oq z 

By the definition of the Wigner function of a pure state described by a wave 
function tp, we have 

W(q,0) = \i;(q)\ 2 

hence 



TT 



1/4 

where is purely real. But, on the other hand, 



Mi) 



tt -i y - 1 1 
W{q,y) = \l aq + bq2 e 4 aq+bqA oc ip(q + -y)ip(q - -y) 
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is only possible if 9 has an imaginary part. 

The equation for the Wigner function of a Yang-Mills gauge field, found by Elze, 
Gyulassy and Vasak, || , uses second quantised quantum mechanics and hence op- 
erate on the mechanical phasespace (q,p) and not on the field theoretic one (^4, n). 
Consequently, they get an infinite order partial differential equation whereas we get 
a finite order functional differential one. We do also, however, get an infinite order 
equation, since for Yang-Mills theory not all equations of motion are constraints, 
only part of them (corresponding to Gauss' law). The remaining equations must 
be found by the same techniques as used by Elze, Gyulassy and Vasak, but this 
time in a purely Hamiltonian framework. 

One should note that for abelian theories, the second term in (|140| ) vanishes and 
we get the classical requirement Q a = (since W ^ in general - it is, a priori, 
possible to have Q a ^ at a point (q,p) provided W(q,p) = of course, hence the 
zeroes of the Wigner function can correspond to points in phasespace at which the 
classical constraints are violated, by continuity of W this can only happen in a 
discrete number of points or in a region disconnected from the rest of phasespace). 
Thus, it is the non-abelianness of the gauged algebra, 5, which leads to quantum 
deformations of the conditions of physical degrees of freedom. 
This corresponds to what was proven in ||] where it was shown that a generalised 
WWM- formalism exists for a large range of algebraic structures, the "deforma- 
tion" of the resulting "classical" phasespace (in particular its curvature), being 
due to the non-commutativity of the generators of the algebra. 
It is not surprising that Yang-Mills theory is anomaly free at this level, since 
anomalies tend to appear through the Dirac operator, and hence through the mat- 



ter fields, [fr^ , 13]. We have only considered Yang-Mills theory in vacuo at this 
stage. The usual anomalies (chiral, parity) should then reappear only when one 
computes the Moyal bracket for the fermion currents. It is already known that the 
WWM symbol of an operator is related to the index of it, and that the Atiyah- 
Singer index theorem normally used to express anomalies is intimately related to 
the entire WWM-scheme, [24]. One should also take notice of the fact that the 
above discussion doesn't take global anomalies into account, such problems are 
beyond the scope of the present paper. 

4.2 Gravitation 

For gravity, we will first consider the ADM constraint algebra, then the Ashtekar 
variables and finally make some brief comments on the newly proposed "Kuchaf 
algebra.", 0, [DJ 
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-i/2 ( 2 _ 7r i ?r j ) + y/ g R= G ijkl ^n kl + JgR (145) 



The ADM constraints, JTT[] , are 

H±(x) = <j 
Hi{x) = 

with gij the 3-metric, it 13 its conjugate momentum, 

{^■(x),7r fc V)}pB = + 5 l «5 l i )5(x,x'), 



(146) 



(147) 



R the curvature scalar of gij (i.e., the three dimensional one) and g the determinant 
of the 3-metric. The first constraint is known as the Hamiltonian one, and the last, 
the TCi, as the diffeomorphism one. The algebra is 



{n ± (x),n ± (x')} PB 

{H ± (x),Hi(x')} PB 
{Hi(x),Hj{x)}-p B 



(g t3 (x)n j (x)+g i3 (x')n j (x'))5 )l (x,x') (148) 
H ± (x)5 A {x,x') (149) 
TLi(x')5j(x,x') +TLj(x)5 ! i(x,x') (150) 



where the subscript, 5i, on the delta function denotes partial derivative with re- 
spect to x 1 . The convention is the standard one in which 5(x, x') is a scalar in the 
first argument and a density in the second (the curved spacetime Dirac 5 has a 
g- 1 / 2 in it). 

The algebra of the diffeomorphism constraint will not be deformed as their form 
is TCi ~ dir + g 2 ir and thus has = 0. The Hamiltonian constraint, however, 
has as well a git 2 as a g 2 n term, and will consequently not have vanishing W3. We 
should thus expect the algebraic relations involving TC± to receive h? corrections 
(but no higher order corrections since no higher powers of ir are present). This is 
precisely what we find. Moreover, the Christoffel symbols and the ^fg contain, in 
a Taylor series, the metric to infinite order, whence we should expect infinite order 
equations to turn up at some stage. 
A straightforward computation yields 



where 



with 



^mnab 



[H ± (x),H ± (x')} M 
k(x,x') = --< 

5 3 H± 
5gij5gki5ir mn 

- G mna b (g tk g 3 



= ih{H ± (x),H ± (x')} PB + h 3 k(x,x' 



(151) 
(152) 

(153) 



2 y 



9nb5 m 5 a + g ma 5 n 5 b - 
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9abSid l n - g mn S% + g bn S l m S k + g am S l n S^ (154) 

= S^6^8g mn (155) 
= -\g mn G ijkl + g- 1 ' 2 { gjl {S?6£ + + to(^ n + to- 

9u8?Sj-9iiSk%) (156) 

One should note that [H±, fcjjy 7^ hence we get an anomaly which is not even a 
central extension of the original algebra. Explicitly 

k] M _ Q PB + «33 hSoi ^_£l_ ^ „ (157) 

For the ADM constraints, the structure coefficients depend on the fields, conse- 
quently the anomaly too depends upon (g, tt). 

Similarly, the relation mixing 7i± and Tii receives a h? correction of the form 

/\ 1, S 3 H ± (x) 5 3 Hi(x>) 

Ki(X, X ) 



8 5iri k 5iT lin 5g ah 5g jk 6gi m 5TT ab 
which one easily finds to be 



with 



5g jk ( X ')5g lm ( X ")5ir ab (y) 



5 c {a 5^(x,y) {6\J%6( X , X ") (-g rk Ti c 6( X , X ')+ 

\o rs {V%dc + 5\ s 5%d r - 5\f r) d s ) «5(*,z')) + 

5U k ) 5(x,x')[-g rm T l rc S(x,x")+ 

+\a rs {s\ r ST)dc + S[ s 6%d r - S\ c S^d s ) - 

<?™<5(x,x") [ 5 ri / m P rc - 5 r V m r< c + 

^/V* (*(A^ C + - 5\ r 5%d s ) - 

\g rl g sm (s{ r 5 k s) d c + 5{f c) d r - 5{ r 5 k c) d s )] 5(x,x')\ + 



(c ^ r,n ^ c,r ^ n) (159) 
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The spatial diffeomorphism subalgebra spanned by Hi does not receive any quan- 
tum corrections since the constraints are only linear in the momentum. 
The set of physical states are defined as the functions W satisfying the two infinite 
order functional differential equations 

= [H±,W]+ (160) 
= [Hi,W]+ (161) 

these are infinite order since the Christoffel symbols (and hence the covariant 
derivative and the curvature scalar) has an inverse metric in them, similarly the 
supermetric Gijki too has an inverse metric inside. Thus the constraints are not 
polynomial in the metric, but instead "meromorphic" . 
Written out more explicitly, the physicality conditions read 



OO i 

= 2H±W + ^(-l) fe 2~ 2fc - 1 /i 2fc ( 
k=i ^ 



5 2k H i 5 2k W 



Sg iljl ...5g i2k _ lhk _ 1 5g mn ...5^-^-^Ti mn 



5 2k H ± 5 2k W 
2jI\j — |— 

Sg iljl ...dg i2k _ lhk _ 1 6TT mn 5^... 5^-^-1 5 g mn 

k(2k _ u ^ ^ W 

Sg hjl ...8g i2k _ 2hk _ 2 5^-^-i8iT™ Stt^i ...5Tr^-^^5g t2k _ lj2k _ 1 5g : 

(162) 

5 2k Hi 5 2k W 



■jnrn 



2HiW + ^2{-l) k 2~ 2k - 1 h 2k ( 



k=1 \ u y«ui •••"y«2fcj2fe 

2k P k H t 5 2k W \ 

Sg iljl ...Sg i2k _ lj2k _ 1 STr^h k S^...S^-^ S g l2kj2k ) 

Since the constraints for gravity in the ADM formalism are non-polynomial the 
equations defining the physical state space become infinite order. If one as- 
sumes the Wigner function to be analytic in h, one can Taylor expand it W = 
H™=o^"^n> an d arrive at the following recursive formulas for the n'th order co- 
efficients, W n 

= H L W = HiW (164) 



i 1 > / / 
= 2H ± W N + ^2{-l) k 2- 2k - 1 ( 
k=l ^ 



6g hjl ...5gi 2k _ l32k _ 1 6ir mn 8^...Si:^-^-^Sg n 



2k S^H± & \ w 

S9nn---Sgi 2k _ 2hk _M 2k - lhk - lS ^ mn ^ tin ---^-^-^g i2k _ lj2k _ 1 5g mn J N 2k 

(165) 
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= 2H i W N +Y(-l) k 2- 2k - 1 - -4-^ = . ■ . . WN ~ 2k . (166) 

where fc, N > 1. These equations are not enough to completely specify the Wigner 
functions. 

The conclusion so far is then that in the ADM formalism gravity is anomalous 
when one attempts a deformation quantisation. The question is, then, whether 
one can lift these anomalies or not. 

Letting TCq = TC± and fi = (0, i) we then want quantum constraints satisfying 

[H^x),H v {x')} M = c^(x,x')H p (x) +dP LV {x,x')H p (x') (167) 

where the structure coeffcients c^,<i^ depend on the coordinates only through 
Dirac delta functions (and their derivatives) and directly through the phasespace 
variables. Because of the complicated nature of the constraints in the ADM for- 
malism I have not been able to find a good set of quantum constraints. 
We saw that the anomalous nature of the quantum deformed algebra of the con- 
traints in the ADM formalism were due to the constraints being non-polynomial. 
It is therefore interesting ot consider another formulation, the Ashtekar variables 



[16 1, where the constraints are polynomials. In this formulation the canonical coor- 
dinates are a complex s«2-connection Af and its momentum (a densitised dreibein) 
E l a , and the constraints are 



H = F«EiEie bc a (168) 
Ga = DiEl (169) 
Vi = F*Ei (170) 

It will turn out that in these variables the anomaly is much simpler, namely merely 
a central extension. 

Since the Ashtekar variables bring out the analogy between general relativity and 
(complexified) Yang-Mills theory due to the isomorphism so(3, 1) ~ SU2 <8> C, we 
can use our knowledge of the deformation quantisation of Yang-Mills systems to 
see that only the following two brackets can receive any quantum corrections, and 
these only to lowest order 

[H(x),H(x')] M = iHn(x)^)hB + ^( J^f E ^f A 



[x <-> x 

(171) 



[H(x),T>i(x')] M = «Mtt(*),Pi(*')}FB - - ^EHA 5 7aHE (172) 
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where we have suppressed the indices on the A, E. An explicit and straightforward 
computation gives 

u 3 (H(x),H(x')) = -12i6(x,x')(Ei(x)A»(x)-Ei(x')A«(x')) (173) 
L0 3 (H(x),V m (x')) = 9iS, m (x,x') (174) 

We notice that the first of these vanish in the sense of distributions, hence the only 
quantum correction is the constant (w.r.t. the phasespace variables) 9i5, m (x, x'). 
Consequently, the anomalous nature of gravity shows itself in the Ashtekar vari- 
ables simply in a central extension of the constraint algebra (similar to a Schwinger 
term in current algebra). 



[H(x),Vi(xf)] M = ih{n(x),Vi(x')}p B - 9ih 3 5i(x,x') 



(175) 



As we have seen earlier such central extensions can be "lifted" by means of a 
redefinition of the quantum constraints involving negative powers of h. 
Furthermore, since the constraints are polynomial in the phasespace variables the 
equations defining the physical state space, C p h ys , become finite order differential 
equations. Explicitly, since the constraints are at most quartic in the phasespace 
variables 



1 



0=[H,W]+ = 2HW--H- 



rpk rpl r bc^a 
a b v a e ( 



S 2 W 
f 5EHE l t 



5 2 W 



5 2 W 



be rpa 

" l HA\bA c 3i 



,5 4 ef 6*W 
+ —n e hr e n J - ? 

4 6E*8E}.8A%8A{ 



0=[Vi,W]+ = 2ViW--h 2 



* a ef EJ a 



5 2 W 
5E l e 5E j f 



5^) + e a mn (5TSfAf + 5^Ar 



= [Ga,W] + = 



2Q a W + -ih 2 5{e c ab — . 

4 k ab 8A c JA k 



f c " 1 f b> 8~E^M 



+ 



5 2 W 
8E k bA a - 



(176) 



(177) 
(178) 



These coupled equations constitute the equations for the Wigner function for 
Ashtekar gravity in vacuum. 

So far the Ashtekar variables have turned out to be much more tractable than the 
ADM-approach; the "anomaly" was simply a Schwinger term and the equations 
for the physical states are finite order. There is one problem, however, which one 
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seldom take into account in this formalism, namely the reality conditions. The 
physical variables have to correspond to a real metric, not a complex one. This 
condition is a second class constraint, |17|j , and this is where the real problems with 
the Ashtekar formalism lies. But since we have seen that second class constraints 
can be treated rather easily in this formalism, the conclusion must be that the 
Ashtekar approach to gravity, all things considered, is the most fruitful one. 
As a final comment, let us digress on the recently proposed "Kuchaf algebra", 



14 1 . Here the Hamiltonian constraint Tl± of the ADM approach is replaced by 
an abelian constraint KL (a scalar density of weight lo), leading to the algebra of 
x Diff(S), i.e., 

{/C(x),/C(x')}pb = (179) 
{K{x),Hi{x')} VB = K,i5{x,x')+u1C5,i{x,x') (180) 

{ni(x),nj(x')}p B = ni(x)5j(x,x')+nj(x')5 )t (x,x') (isi) 

There is an appealing interpretation of this algebra in terms of fibrebundles over a 
three-manifold, pq| , but since the only known representations of this algebra are 
the ones formed from the ADM constraints, such as 

K, = TCj_ — JiiHjQ 1 - 1 

and even more non-linear expressions (the general solution being found in |lj 
we will not be able to compute the Moyal brackets in any satisfactory way - they 
will a priori involve infinite order differentials. There is, however, still the hope 
that some of the solutions, i.e., some of the explicit expressions for the /C's, will 
turn out to simplify the Moyal bracket and will hence be somehow favoured. At 
the present, it has to be admitted, though, that this is a rather faint hope. 
It is not enough for a constraint algebra to be nice, it is also important for the rep- 
resentation of the constraint algebra as functions on phasespace to be sufficiently 
simple, in order for the deformation quantisation scheme to be really tractable. 
This of course holds in practice for any quantisation scheme, but contrary to some 
other schemes, deformation quantisation is at least in principle able to handle 
constraints of arbitrary complecity - it is merely a matter of computational con- 
venience (or laziness). 

In a sequel paper we will concentrate on the deformation quantisation of gravity, 
and we will interpret the set of quantum constraints in the Ashtekar variables, 
find the relationship to the loop formalism and knot invariants and finally find an 
explicit solution related to topological field theory, [28|. 
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5 Underlying Geometrical Structure 



In the present picture one quantizes a classical theory, as described by an algebra 
of observables .Ao = C°°(T), by deforming it to obtain a noncommutative algebra 
An ~ Ao <8> C((h)) (isomorphism as vectorspaces, not as algebras), with a new, 
twisted, product / ', g i— > / * g = f g + 0(h) . 

We can consider An as a "field" (in the language of C*-algebras) over the real 
axis, parametrised by h G R. This is the point of view taken by Nest et al, 24]. 



Alternatively, we can consider a sheaf, [g9|, |3Cj , A over the topological space X 
where the stalks are given by 



A a 



x < 
A x = 
An x > 



(182) 



where h is taken to be a free parameter, h := x. This is obviously a fine sheaf, 
since a partition of the identity trivially exists, furthermore, it is a sheaf which is 
almost constant. This expresses quantisation as a sheafication process. 
We have already commented a bit on some underlying cohomological structure, 
here we will elaborate on this. First, each / € Ah can be written as a Laurent 
series (we can without loss of generality assume only one negative power of h to 
be present) 

f = f-ih- 1 + h + hh + f 2 h 2 + ... fieAo (183) 

This leads to the definition of a series of natural differentials, each with their own 
interpretation 



5 


f- 


+ /o 


s q 


f- 


-hh + h 2 f 2 + 


5- 


f- 




5+ 


f- 





where 5± are inspired by the corresponding definition for the Weil complex, |2J 
Clearly 5q = S 2 , = = so these three all define cohomology theories by consid- 
ering the trivial complex 

O^Ah^Ah^Ah^- (184) 
where 5 denotes any of the differentials S±,5q. The cohomologies are trivial 

H°(8 ) = Ao H°(8±) = 5±A h (185) 
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thus, H°(5q) gives the classical algebra back, whereas H°(5—) gives the anomalous 
part. 

A more interesting complex is 

O^Ah^Ak^Ak^... (186) 

with alternating So, 5 q . This is indeed a complex since obviously So5 q + 5 q 5o = 0. 
This has the same cohomology as the 5o- com pl ex i since Im<5o = Ker<5 g . 
But, instead of considering the cohomology of a single stalk, it is better to consider 
the cohomology of the entire sheaf. This will also allow us to take the physicality 
condition into account. This condition can be rewritten as 

S / := [rf <W] M = (187) 

where now r] a a constant Grassmann numbers and the Moyal bracket is graded 
appropriatly - when / is Grassmann even this becomes the anti- Moyal bracket. 
Then, s 2 = 0, and s is very similar to the BRST differential. By extending the 
algebras A x to include the constant Grassmann parameters r/ a , we get a complex 
of sheaves 

s : A q -> A q+1 (188) 
where A q denotes the sheaf obtained from A by tensoring 

A q = A® G®G® ...®& (189) 
s 

where G denotes the space of the ^-parameters. This leads us naturally to the 



concept of hyper cohomology, |30|]. Let {J} be a covering of X = M by open sets, 
i.e., {/} consists of open interval ]a,b[ such that their union gives the entire real 
axis. Given such a covering we can define Cech cochains with values in the sheaf 
complex .A*, the set of these is denoted by C P (I,A*) and is defined by, p9j] , 

C P (I, A*) = {/ : I p+1 -> A*} (190) 

The Cech cochains come with a differential, 8, given by 

6f(i , ...,I P+1 ) = /(J n l p+1 ,h, ...,l p ) + ^(-i)*+i/(/ , .... 4 .... 4 ...i, n ij) 

i>j 

(191) 

We hence have a bicomplex, C p (I,A q ), the cohomology of the total complex is 
denoted by simply H*(I,A*). Now, any refinement of the cover induces homo- 
morphisms among the cochains and hence also among the cohomology modules, 
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consequently, we can perform the limit of ever finer coverings to obtain the hyper- 
cohomology 

M*(X,A*) := limH*(I,A*) (192) 
Since the sheaf is fine, we have 

U p {X,A q ) = p^O (193) 

On the other hand 

B°(X,A q ) = limKer(<5 + s) (194) 

In ghost number zero, we therefore have 

U°(X,A°) =C phys (195) 

This method is compatible with BRST-techniques. One simply defines the quan- 
tisation by means of the abovementioned sheafication procedure. This remedies 
the usual shortcoming of canonical quantisation techniques, such as the BRST, by 
avoiding finding an operator realisation. One just considers the BRST symmetry 
as a classical symmetry, and then deforms it by replacing the BRST-complex by 
its sheafication over the real axis. 

6 Conclusion 

We have studied the deformation (i.e., Moyal) quantisation of Hamiltonian systems 
with constraints. Instead of looking for an operator correspondence qi,p l — > qi,p l 
with the usual rule {/, g}pB — ► ^[/>ff]) we keep the classical phasespace but replace 
the Poisson bracket with the Moyal bracket, \f,g\M = 2^{/><?}pb + 0(h 2 ). It is 
know that this is always possible. By introducing sufficiently many phasespace 
variables and constraints one can assume the classical phasespace to be flat (i.e., 
to have a global patch of Darboux coordinates, {qi,p^}pB = Sj). 
A priori one have to replace the classical constraints (j) a by quantum deformed 
versions <3? a = ^ n h n &a^ with = (j) a in order to keep the constraint algebra 
when one replaces Poisson brackets by Moyal ones. We saw that for Yang-Mills 
systems we could take $ a = <j> a , but for systems where the constraints are more 
than cubic in the phasespace variables, one would in general have ^ (pa- 
Second class constraints could be handled by a similar device, turning them into a 
an Abelian algebra under Moyal brackets. This was illustrated by the case of the 
O(N) non-linear cr-model. 

The dynamics of a constrained Hamiltonian system was captured by the Wigner 
function, where the classical condition <p a = was replaced not by the operator 
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equation <p a \^) = but by the algebraic condition [<I> a ,VF]^ = resembling the 
BRST approach. 

Various problems can arrise: (1) the appearance of second class constraints, (2) the 
appearance of first class constraints not in involution with the Hamiltonian, (3) the 
appearance of anomalies. It has been shown in this paper how all these problems 
can be remedied in deformation quantisation, provided one is willing to include 
negative powers of the deformation paramater h in the expansions of the quantum 
objects and to let the structure coefficients receive ^-corrections. Hence, upon 
replacing An = A®C[[h]] („4- valued formal power series in h) by An = A®C((h)) 
(the set of A- valued formal Laurent series in h). This of course implies that the 
naive classical limit h — ► becomes singular. The correct classical limit appears 
only upon taking principal values. Thus, the naive classical limit A c i as = lim^o A 
is replaced by 

A clas = PP lim A = 5 A (196) 

h— >o 

where PP denotes principal part and where 5q is the differential defined in the 
previous section. 

With this definition the classical part is still the frP component, there is just a 
subtle difference between picking out this component and taking the limit h — > 0. 
With this caveat, deformation quantisation emerges as a very powerful quanti- 
sation scheme indeed. In particular, on should note the very explicit form the 
relations determining the ^-modification to the classical objects take - an explicit 
form allowing us in certain examples to arrive at an explicit solution to all orders 
in H. Consequently, deformation quantisation is a very constructive approach to 
quantisation. The caveat about the principal part also illuminates the nature of 
the classical limit, and in particular how ill-behaved classical theories can arrise 
from well-behaved quantum ones due to the singularity of the limit. Alternatively, 
one can view this as illustrating the quantum "smearing out" of problems in the 
classical theory. With this in mind, it would be very interesting to study the ap- 
plication of deformation quantisation to quantum chaos. 

Yang-Mills theory was then quantised in this scheme and we saw that in the ab- 
sence of matter, the theory was anomaly free. 

Gravitation was treated in three different manners, first the standard ADM ap- 
proach, secondly in the Ashtekar variables and finally with a recently proposed 
new constraint algebra. We found that gravity in the ADM formalism was anoma- 
lous and lead to infinite order equations for the physical states - both problems 
stemming from the non-polynomial nature of the constraints. The Ashtekar vari- 
ables, however, turned out to aquire merely a central extension, which can be 
lifted. Moreover, since these constraints were polynomial, the equations picking 
out physical states became finite order. 
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